Persistence in systems with algebraic interaction 
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Persistence in coarsening ID spin systems with a power law interaction r _1_<T is considered. 
Numerical studies indicate that for sufficiently large values of the interaction exponent a (a > 1/2 
in our simulations), persistence decays as an algebraic function of the length scale L, P(L) ~ L~ . 
The Persistence exponent 6 is found to be independent on the force exponent a and close to its 
value for the extremal (a — > oo) model, 9 = 0.17507588 .... For smaller values of the force exponent 
(a < 1/2), finite size effects prevent the system from reaching the asymptotic regime. Scaling 
arguments suggest that in order to avoid significant boundary effects for small a, the system size 
should grow as [O (1 / a)] 1 1 17 . 
PACS numbers: 02.50.Ga, 05.70.Ln, 05.40. +j 



Coarsening dynamics of one-dimensional systems with 
a power-law V(r) ~ r~ a ~ 1 interaction between spins has 
recently been studied by Lee and Cardy [Q , and Ruten- 
berg and Bray It had been established that after 
quenching from a high-temperature disordered phase to 
T = these systems develop a domain structure charac- 
terized by a single lengthscale L(t). A naive argument 
based on the law of motion for domain walls, L ~ 
(where L~ a is a typical force between domain walls), pro- 
duces an asymptotically correct time dependence of L, 



L(t) ~ i 1/1+<T . 



(1) 



Other properties of this system, including correlation 
functions and domain size distribution, have been studied 
m |] as well. 

In this paper we shall look at another facet of ID phase- 
ordering systems with a power-law interaction: what 
fraction P of spins have never changed sign up to the 
time t? Or, equivalently, what fraction of the space has 
never been crossed by a domain wall? Such a property of 
coarsening systems is usually called persistence and has 
recently become a major subject of research in statistical 
physics • Let us briefly review some known results 
in this field relevant to our problem. In || the exact 
solution was found for persistence in an ordering system 
with extremal dynamics. The extremal dynamics limit 
is reached by formally setting a — > oo, which means that 
interactions become infinitely short-range. In this limit, 
coarsening proceeds by consecutive shrinking and disap- 
pearance of the current smallest domains in the system, 
while other domain boundaries remain virtually motion- 
less. It was established in || that persistence at a stage 
of evolution when the average domain size is L is propor- 
tional to L e , where the exponent — 0.17507. . . is the 
solution of an implicit integral equation. 

In H persistence exponents have been calculated for 
coarsening ID Potts models with Glauber dynamics. For 
the 2-state Potts (Ising) model, persistence decays as 
t~ e , — 3/8, or in terms of the average domain size 
L, P(L) ~ L- 3 / 4 . 



The following conclusion can be drawn from a compar- 
ison of persistence exponents for extremal and Glauber 
dynamics. Extremal dynamics is more efficient in pre- 
serving persistence, since the motion of domain walls is 
always directed towards their ultimate annihilation part- 
ners, while in the case of Glauber dynamics, domain 
walls perform random walks and sweep through a larger 
amount of space, which otherwise could have remained 
persistent. The extremal dynamics exponent 6 sets a 
lower bound on persistence exponents for systems with a 
finite force exponent a. It is easy to visualize a scenario 
when a domain wall first moves away from its ultimate 
annihilation partner, and then, after the stronger force 
source disappears, it turns back. Such events result in 
spin flips on parts of the line that belong to a surviv- 
ing domain and would have been left untouched in the 
extremal dynamics case. 

The results presented below suggest that this lower 
boundary 9 = 0.17507 ... is in fact the exact value of the 
persistent exponent for arbitrary a > 0. 

Let us formally introduce our model: We consider 
coarsening of the ID 2-state spin system with a long- 
range ferromagnetic Hamiltonian: 



H 



a 4>j (xi- Xj) 



cr+l ' 



(2) 



After quenching from a high-temperature random phase 
to T = 0, coarsening dynamics for this system is deter- 
mined by the motion of domain walls, governed by the 
Langevin equation. The velocity of a wall is equal to 
the sum of pairwise forces from other walls, with walls of 
the same signs repelling and walls of the opposite signs 
attracting each other: 



~dt 



F 



i+J 



sign(ri - rj)F i:j 



n 



(3) 



(4) 



1 



When the adjacent walls meet, they annihilate. As we 
mentioned, the degree of coarsening is uniquely charac- 
terized by a typical domain size L(t) ~ i 1 /( 1 + tr ). We 
measure the fraction of space P(L) that has never been 
crossed by a single domain wall as a function of this 
lengthscale L(t). We perform molecular dynamics simu- 
lations of the model for the following values of the force 
exponent: a = 3/2, 5/4, 1, 3/4, 1/2, 1/4. Each run starts 
with a system consisting of No = 100000 domain walls; 
results for each a are averaged over 20 initial configura- 
tions. Open boundary conditions with no replicas added 
to the boundaries are used. To speed up the evaluation 
of forces, a ID multipole expansion has been performed, 
and terms of up to quadrupole order were taken into ac- 
count ||. 



and get inside the territory of the future survivor, causing 
some excessive spin flips. 
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Fig. 1. The log-log plot of persistence P(L) vs. average do- 
main size L for various force exponents a. The straight line 
corresponds to P(L) ~ L~ e . 

The results for persistence as a function of the average 
domain length L are presented in log-log form in Fig. 1. 
Except for small force exponents (a — 1/4 and later evo- 
lution stages for a = 1/2) all the curves collapse at a 
line with a slope ~ —0.175, which corresponds to a = oo 
extremal model. 

Statistical error bars are shown in Fig. 2 for a single 
set of data (cr = 5/4). 

Our simulations suggests that scaling of persistence, 
corresponding to a = oo, is valid for all other not very 
small a. The following asymptotic argument helps to un- 
derstand why this is so. At any current moment of time, 
persistent spins are mostly contained in the domains that 
were expanding at almost all previous stages of coarsen- 
ing; i.e. these domains were larger than the average at 
those stages. If one of such large domains is surrounded 
by two small neighbors, it would most probably grow out- 
wards, and no spin flips, additional to those inevitably 
caused by directed coarsening itself, would happen. The 
situation may be different if two or three big domains are 
adjacent to each other: their domain walls may wonder 



-1.0 1 1 1 1 1 1 1 1 1 1 1 

0.0 1.0 2.0 3.0 4.0 5.0 

lg(L) 

Fig. 2. The log-log plot of persistence P(L) vs. average do- 
main size L force exponents a = 5/4 with statistical error 
bars. The straight line corresponds to P(L) ~ L~ e . 

We can estimate the characteristic scale of such a 
persistence- loosing event. A typical distance AL that 
a wall of domain of size Lq, surrounded by a group of 
domains of similar sizes, travels during time t is 

AL ~ L - t?' x+ ° « L(t)[^r. (5) 

La 

Here L(t) ~ t 1 / 1+cr is the average domain size at time 
t. For positive cr, AL becomes negligible compared to 
L when 3> 1, hence the number of spin flips addi- 
tional to those present in extremal dynamics coarsening 
scenario becomes negligible. Another conclusion that fol- 
lows from Eq. (||) is that for small a, the crossover time to 
P(L) ~ L 6 scaling must be larger since the system must 
develop structure that includes sufficiently large domains. 

However, besides long intitial transitional times, there 
is another reason for the breakdown of scaling for small cr 
that we observed in our simulations. Let us first consider 
the opposite to a — oo case of a = 0. In this limit forces, 
are distance- independent, and domain wall dynamics (|^) 
is described by the equation 

5=X)(-l) <+J 'sign(r i -r J ). (6) 

Since domain walls come in pairs, the sum in Eq. (^) is 
equal to ±1. That means that all walls have the same 
constant velocity with odd-number walls moving to the 
left and even- number moving to the right. The whole 
system becomes a collection of independently collaps- 
ing and growing domains. This clearly violates the scal- 
ing (Q); in fact the a = system has two lengthscales, 
Lq — 2vt and Lq + 2vt, where Lq is the average initial 
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domain length and v = 1 is the velocity of domain walls. 
For an exponential distribution of initial domain sizes, 
W(Lq) = exp(— Lq), persistence can be expressed as 



P(t) 



J °° W{x)dx + /~ ( W(x)dx _ 1 + exp[-2t] 
2 J °° W(x)dx ~ 2 



(7) 



Systems with few particles and small a coarsen almost 
according to the a — scenario: particles across the 
whole system feel the presence of the boundary. Odd 
and even -number walls tend to move predominantly to 
the left and right, respectively, independent of the posi- 
tion of their nearest neighbors. 

To probe whether the deviation from P(L) ~ L e scal- 
ing in persistence behavior is caused by a — finite size 
effects, we do the following measurements. First, for the 
system of the same initial size (N — 10 5 ) we plot the 
average domain length L(t) as a function of time and 
compare it to the L ~ t 1 / 1+a prediction. 



that our system is never in the true infinite-size regime 
for a = 1/4, finite size effects are becoming evident for 
a = 1/2, even at early stages of evolution, and only for 
a > 3/4 the boundary effects could be neglected. 
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Fig. 3. The log-log plot of average domain size L(t) vs. time 
for various force exponents a. Straight lines correspond to 
scaling predictions, L(t) ~ t 1 ' 1+a . 

Results for this simulation are presented in Fig. 3. One 
can see that the system with a = 1/4 is never in scaling 
regime ([!]), system with a = 1/2 behaves according to 
(|l|) only up to some intermediate stage of evolution. For 
all other force exponents a > 1/2, for a certain period 
of evolution after short transitional time, typical domain 
sizes scale according to (0). Another check of whether 
a system feels the presence of the boundaries and there- 
fore crosses over to the a = coarsening regime, is to 
measure directly the fraction of domain walls B(L) that 
move opposite to the direction prescribed by boundary 
effects. In Fig. 4 we plot the fraction of even-number do- 
main walls moving to the right and odd-number domain 
walls moving to the left for the systems initially consist- 
ing of the same number of domains, N = 10 5 . For finite 
a > and a truly infinite system this fraction should 
be equal to 1/2, for a = it should be 0. We observe 



Fig. 4. Number of domain walls B(L) that move opposite 
to the direction prescribed by boundary effects vs. average 
domain length L for various force exponents a. 

Finally, we present a rough estimate of how big a sys- 
tem should be for a particular value of a <C 1 to avoid 
significant finite-size effects. We evaluate a force F1-2, 
exerted on a test domain wall by a dipole pair of neigh- 
boring domain walls, 



F1-2 



(j)" -(^)° ^(|)Vm2 + 0(a 2 )) (8) 



and compare it to the force Fn exerted on the same test 
domain wall by a single domain wall at the the system 
boundary. 



F N 



2 c 



O) 



Here L and ./V are the typical domain length and the 
number of domains in the system. The boundary effects 
become significant when these forces are of the same or- 
der. Hence for particular a, the minimum number of 
particles to avoid finite size effects N m i n is 



2( 



1 



cr ln2' 



(10) 



We have observed (see Fig. 4) that for a = 1/2, N min w 
10 5 . Assuming the following parameterization of minimal 
size of the system for small a, N m i n — (constant/cr) 1 ^ 
fitting it to N min {<7 = 1/2) = 10 5 , we obtain 
,,(<t = 1/4) w 1.6 x 10 u . This is well beyond the 



and 



limits of computational power available to us. 

In summary, we presented numerical evidence and a 
scaling argument suggesting the universality of persis- 
tent exponent for extremal model, 9 = 0.17507588..., 
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for models with arbitrary force exponents a > 0. We 
found that a deviation from scaling for persistence, that 
happens for small a, is accompanied by a similar de- 
viation from scaling for a typical domain size L{t) and 
is caused by finite size effects that cause crossover to a 
(7 = coarsening scenario. We estimated that in order 
to avoid boundary effects, the system size should grow as 
[C(l/cr)] lcr . A possible extension of this work is for higher 
dimensional systems, though the duality between domain 
walls and spin dynamics that was extensively used for this 
work, may not be so straightforward to apply. 

The author would like to thank P. Krapivsky, A. 
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